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SUMMARY 

Ames Research Center has the lead role 
among NASA centers to conduct research in 
computational fluid dynamics. The past, the 
present, and the future prospects in this 
field are reviewed. Past accomplishments 
include pioneering computer simulations of 
fluid dynamics problems that have made com- 
puters valuable in complementing wind tunnels 
for aerodynamics research. The present 
facilities include the most powerful com- 
puters built in the United States. Three 
examples of viscous flow simulations are 
presented: an afterbody with an exhaust 

plume, a blunt fin mounted on a flat plate, 
and the Space Shuttle. The future prospects 
include implementation of the Numerical Aero- 
dynamic Simulation Processing System that 
will provide the capability for solving the 
viscous flow field around an aircraft in a 
matter of minutes. 

INTRODUCTION 

A major goal of the National Aeronautics 
and Space Administration (NASA) is to conduct 
an effective and productive aeronautics 
research and technology program that contrib- 
utes materially to civil and military avia- 
tion. A key objective in attaining this goal 
is the advancement of both the understanding 
of and the capability to predict aerodynamic 
phenomena through the use of wind tunnels and 
high-speed computers. The results of past 
efforts are exemplified by the large number 
of commercial airliners built in the United 
States and operated worldwide, including 
Japan (fig. 1). 

In support of NASA's goal, Ames Research 
Center has established and strives to main- 
tain one of the world's premier aeronautical 
research establishments. In computational 
fluid dynamics, Ames Research Center has the 
lead role among NASA centers and has 



Figure 1.- Japanese commercial airliners in 
1983 (number in service). 


installed the most powerful computers built 
in the United States. The research program 
encompasses the development of numerical 
methods and computer codes for fluid-flow 
phenomena, experiments for verifying the 
codes, and application of the codes to prac- 
tical aerodynamics problems. This paper 
reviews the past, the present, and the future 
prospects of computational fluid dynamics, 
based on experiences at Ames Research Center. 

THE PAST 


The history of computational fluid 
dynamics at Ames Research Center began 
25 years ago with the establishment of NASA 
when Ames Aeronautical Laboratory was trans- 
ferred from the National Advisory Committee 
for Aeronautics (NACA). Van Dyke (ref. 1) 
used an IBM 650 computer to obtain numerical 
solutions of the inviscid flow field around a 
blunt body traveling at supersonic speeds. 
Reentry into Earth's atmosphere by space 
capsules and missiles was the application of 
such studies. This era also saw the begin- 
ning of the commercial jet age, with trans- 
continental and transatlantic flights by 
Boeing 707 and Douglas DC-8 airliners. 

Fifteen years ago, computational fluid 
dynamics became recognized as a separate 
field of study. The combination of the 
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THE PRESENT 


availability of the IBM 360 series computers 
and development of new algorithms allowed 
solution of problems previously considered 
Infeasible. MacCormack (ref. 2) devised an 
algorithm for the Navler-Stokes equations and 
solved numerically the interaction between a 
shock wave and a laminar boundary layer; 
several investigators obtained solutions of 
the small-disturbance equations for inviscid 
transonic flow over an airfoil. These 
results were reported at the Second Interna- 
tional Conference on Numerical Methods in 
Fluid Dynamics held at Berkeley, California, 
September 15-19, 1970 (ref. 3). Air travel 
was augmented by the introduction of the 
jumbo Boeing 747, now the standard airliner 
for intercontinental flights. 

Ten years ago, transonic flow solutions 
were obtained for inviscid flow over wings by 
Lomax et al. (ref. 4) and for viscous flow 
over airfoils by Deiwert (ref. 5). The Con- 
trol Data Corp. (CDC) 7600 became available 
for fluid dynamics simulations, and the 
ILLIAC IV, a parallel processing computer 
with 64 separate processing elements, was 
installed at Ames. The ILLIAC IV computer 
was used extensively to study turbulence from 
first principles, beginning with the direct 
numerical simulations by Rogallo (ref. 6). A 
major event in air travel was the introduc- 
tion of the British-French Concorde supei — 
sonic airliner. 

Five years ago, unsteady viscous flow 
solutions were obtained for the aileron buzz 
problem by Steger and Bailey (ref. 7), and 
steady three-dimensional viscous flow solu- 
tions were obtained for axisymmetric bodies 
at angle of attack by Pulliam and Steger 
(ref. 8) and by Hung (ref. 9). In a program 
conducted by Stanford University and Ames 
Research Center, calculations of turbulence 
were performed on the ILLIAC IV computer, 
using the large-eddy simulation technique, in 
which large-scale eddies are calculated from 
the unsteady Navier-Stokes equations and 
subgrid-scale eddies are modeled. Results 
for turbulent channel flow were obtained by 
Moin and Kim (ref. 10). During the past 
5 years, the solution of three-dimensional 
viscous flow problems has become the major 
focus of interest. 


The status of computational fluid dynam- 
ics at Ames Research Center is that of a 
maturing science that complements the role of 
wind tunnels in aerodynamics research. The 
computer facilities will be described first, 
followed by a short discussion of numerical 
methods and by the presentation of three 
examples of viscous flow simulations: an 

afterbody with an exhaust plume, a blunt fin 
mounted on a flat plate, and the Space 
Shuttle . 

Computer Facilities 

The computers used primarily for compu- 
tational fluid dynamics research at Ames 
Research Center include Cray X-MP/22 and CDC 
Cyber 205 supercomputers and several 
VAX-1 1/780 minicomputers. All the computers 
are linked by high-speed lines so that data 
can be transferred between any nodes in the 
network as shown in figure 2. 



Figure 2.- Ames Research Center 
supercomputer network. 


The Cray X-MP supercomputer has two 
processors and storage totaling 2 million 
64-bit words of main memory, 16 million words 
of extended solid-state memory, and 1.2 bil- 
lion words on disks. The Cray is operated in 
a batch mode, with jobs submitted through the 
Cyber 720 front-end processor and through 
four VAX systems connected by Hyperchannel 
links. The current Cray operating system 
does not allow multitasking, but a future 
version will; multitasking allows a job to be 
executed on both processors. 
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The Cyber 205 supercomputer has four 
pipelines and storage totaling 4 million 
64-bit words of main memory and 1.2 billion 
words on disk. This computer is also oper- 
ated in a batch mode, with jobs submitted 
through the Cyber 835 front-end processor and 
through two VAX systems connected through a 
loosely coupled network. 

A typical VAX-1 1/780 system connected 
directly to the supercomputers has 8 mega- 
bytes of main memory and in excess of 
600 megabytes of disk storage. Individual 
researchers work interactively with the VAX 
system from a cathode-ray-tube terminal to 
edit programs and submit jobs to either the 
Cray X-MP/22 or Cyber 205 supercomputer. The 
results are processed on the VAX system and 
may be displayed on the cathode-ray-tube 
terminal, as shown in figure 3. 



Figure 3-~ Cathode-ray-tube terminal used to 
access Ames computer facilities. 

Numerical Methods 

The levels of approximation for the 
governing equations solved in current compu- 
tational fluid dynamics studies are shown in 
figure 4. Early work solved the inviscid 
equations, beginning with linearized theory, 
followed by the small-disturbance and full- 
potential equations, and, most recently, the 
Euler equations. Viscous effects are solved 
separately using the boundary-layer equa- 
tions, which may or may not be coupled with 
the inviscid equations. This approach is 
valid for high-Reynolds-number flow with 
little or no flow separation. 

Recent work has concentrated on the 



Figure 4.- Levels of approximation for fluid 
dynamics equations. 

Reynolds-averaged Navier-Stokes equations, 

with the parabolized approximation allowing 
marching in the streamwise direction for 
supersonic flow and the thin-layer approxi- 
mation neglecting viscous diffusion along the 
surface. Turbulence modeling is required in 
all these studies, with the Baldwin-Lomax 
model (ref. 11) being the one most commonly 
used. The large-eddy simulation method 
requires only that the subgrid-scale turbu- 
lence be modeled. The ultimate goal is to 
solve the complete, unsteady Navier-Stokes 
equations without any turbulence modeling for 
complex configurations. In general, however, 
as the level of approximation rises, the 
geometry must be less complex because of 
limitations on computer speed and storage. 

The selection of the grid is an impor- 
tant step. Some of the criteria to be met 
are that the grid should conform to the body, 
be clustered where flow gradients are steep, 
match the flow physics, simplify application 
of boundary conditions, and produce well- 
ordered equations adaptable to vector pro- 
cessing on the latest supercomputers . The 
grid-generation methods require solution of 
either algebraic or partial differential 
equations; both are in common use. 

Originally, the most popular algorithm 
for the solution of the Navier-Stokes equa- 
tions was the MacCormack method (ref. 2), 
which was explicit in time and used central 
differencing in space. Considerable computer 
time was required, especially for hlgh- 
Reynolds-number flows, because the stability 
condition for explicit methods allows only a 
small time-step for a small mesh. Implicit 
methods such as the Beam-Warming method 
(ref. 12), which allow much larger time- 
steps, were developed as a result and are 
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widely used today. Implicit methods, how- 
ever, require more computer time per time- 
step and generally are more difficult to 
program. Recently, MacCormack (ref. 13) 
modified his original method to make it 
implicit and more efficient by several orders 
of magnitude. 

Example 1: Afterbody with an Exhaust Plume 

The first example of a viscous flow 
simulation Is for an afterbody with an 
exhaust plane. Significant drag can result 
from the interaction between the exhaust and 
the flow field surrounding a missile, rocket, 
aircraft afterbody, or nacelle. Solutions of 
the thin-layer form of the Navler-Stokes 
equations were obtained by Deiwert (ref. 14) 
for flow over an axisymmetric afterbody con- 
taining a centered propulsive jet. Nakahashi 
and Deiwert (ref. 15) have recently improved 
the efficiency and accuracy of the solution 
method by introducing an adaptive grid, which 
clusters grid points in the vicinity of shock 
waves and slip lines where the density gra- 
dients are large. The original grid and the 
final solution-adaptive grid are shown in 
figure 5(a) for a free-stream Mach number of 
2.0, a jet exit Mach number of 2.5, a 
Reynolds number based on body diameter of 
0.65x10®, and an exhaust jet to free-stream 
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Figure 5.- Afterbody flow field with exhaust 


static pressure ratio of 6.0. The density 
contours obtained with the fixed original 
grid and the adaptive grid are shown in 
figure 5(b). Significant improvements in the 
solution are noted in the vicinity of the 
compression shock, barrel shock, and Mach 
disk . 

Example 2: Blunt Fin Mounted on a Flat 

Plate 

The second example of a viscous flow 
simulation is that of a blunt fin mounted on 
a flat plate. Loss of control effectiveness 
and severe local heating are some of the 
consequences of supersonic flow over bodies 
in regions where wings or fins are 
attached. The bow shock wave ahead of the 
blunt fin interacts with the boundary layer 
on the flat plate (or body) and causes sepa- 
ration of the flow as shown in figure 6. 
Solutions of the thin-layer form of the 
Navier-Stokes equations were obtained by Hung 
and Buning (ref. 16) for a free-stream Mach 
number of 2.95, a unit Reynolds number of 
63x10®/m, and a fin diameter of 1.27 cm. 
Particle paths are shown in figure 7 near the 
plane of symmetry for the separated region. 
Particles near the flat plate separate from 
the surface, reverse direction, and spiral 
out past the blunt fin. Particles away from 
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plume. (a) Grid; (b) density contours. 
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Figure 6.- Blunt fin mounted on flat plate. 



Figure 7.- Particle paths near plane of 
symmetry of separated region ahead of 
blunt fin. 


the flat plate are deflected around the blunt 
fin. 

Example 3: Space Shuttle 

The last example of a viscous flow simu- 
lation is for the Space Shuttle. During 
reentry into the atmosphere at high angle of 
attack, intense heating occurs on both the 
windward surface and regions of the leeward 
surface where flow reattaches after separa- 
tion. Solutions of the parabolized Navier- 
Stokes equations were obtained by Chaussee 
et al. (ref. 17) for the flow field around 
the forward two-thirds of the Space Shut- 
tle. Simulated oil-flow lines for turbulent 
flow on the surface are shown in figure 8 for 
the following wind-tunnel conditions: a 

free-stream Mach number of 7.9, a Reynolds 
number of 1.37x10® based on vehicle length, 
and an angle of attack of 25°. Separation is 
indicated on the strake-wing by the coales- 
cence of oil-flow lines, followed by reat- 
tachment part way up the fuselage, and sep- 
aration again on the upper surface. 



Figure 8.- Simulated oil flow lines on Space 
Shuttle surface. 


THE FUTURE 

To ensure continuing leadership in com- 
putational fluid dynamics and related disci- 
plines, NASA has initiated the Numerical 
Aerodynamic Simulation (NAS) Program to 
establish a national computational facility 
at Ames Research Center. The NAS Processing 
System will comprise six subsystems: high- 

speed processor, mass storage, support pro- 
cessing, long-haul communication, graphics, 
and work station. All six are connected by 
the data network subsystem, as shown in 
figure 9. 

The heart of the system will be the 
high-speed processor subsystem, which will 
employ the fastest processors available at 
any given time. The initial configuration, 
which will be operational in 1986, will have 
a Cray-2, with a sustained speed of 250 mil- 
lion floating-point operations per second 
(mflops) , a main memory of 64 million 64-bit 
words, and 15.6 gigabytes of disk storage. 

In 1988, this processor will be supplemented 
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Figure 9.- Numerical aerodynamic simulation 
processing system. 









by one expected to have a sustained speed of 
1 billion floating-point operations per sec- 
ond (gflops) and a main memory of 256 million 
words. The mass storage subsystem will pro- 
vide on-line and archival storage with an 
initial capacity of 200 gigabytes. The sup- 
port processing subsystem will provide inter- 
active support for remote users and 
input/output devices, such as a printer, film 
recorder, and tape drives. The long-haul 
communication subsystem will provide access 
to the NAS Processing System from remote 
sites, Including other government labora- 
tories, universities, and industry. The 
graphics subsystem will provide hard-copy 
graphical output for all users and dynamic 
high-resolution color displays for on-site 
users. The work-station subsystem will pro- 
vide access to the NAS Processing System for 
on-site research scientists and engineers. 

The initial configuration will have 25 self- 
contained units, each equipped with a pro- 
cessor, color graphics display, dot matrix 
printer, and 500 megabytes of disk storage to 
develop computer codes, process output data 
from the high-speed processor, and prepare 
reports . 

Successful implementation of the NAS 
Processing System in 1988 will provide the 
capability of solving the Reynolds-averaged 
Navier-Stokes equations for an aircraft con- 
figuration in a matter of minutes. 

CONCLUDING REMARKS 

Continuing progress in computational 
fluid dynamics will require research and 
development in turbulence modeling, numerical 
methods, and computer hardware. 

The large range of length scales associ- 
ated with turbulent flow precludes solution 
of the unsteady Navier-Stokes equations for 
complex configurations with the computers 
that will be available in the near future. 
Turbulence modeling will be required to solve 
the Reynolds-averaged Navier-Stokes equations 
for three-dimensional flows with separa- 
tion. Both wind-tunnel and computer experi- 
ments are needed. New experimental tech- 
niques using nonintrusive laser devices will 


be complemented by computer simulations of 

turbulence. 

Numerical methods have become more effi- 
cient during the pa3t 15 years, but more 
gains are possible. The need for hundreds 
and, sometimes, as many as thousands of iter- 
ations to solve a steady flow problem is not 
acceptable. The use of implicit methods will 
avoid the time-step restrictions for stabil- 
ity inherent in explicit methods. 

Advances in computer architecture and 
semiconductor technology are resulting in 
supercomputers that are faster and that have 
larger memories (fig. 10). Increasing the 
number of processors is a straightforward way 
to increase the speed, and the processors may 
be assigned different tasks for the same 
problem. The numerical method must take into 
consideration the computer architecture in 
order to optimize the use of the available 
hardware . 



Figure 10.- Peak processing speeds for 
current and future supercomputers. 
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